3.870 \(\int \frac{x^2}{(a+b x^2+c x^4)^2} \, dx\)

Optimal. Leaf size=221 \[ -\frac{x \left (b+2 c x^2\right )}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}+\frac{\sqrt{c} \left (2 b-\sqrt{b^2-4 a c}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (\sqrt{b^2-4 a c}+2 b\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{\sqrt{b^2-4 a c}+b}} \]

[Out]

-(x*(b + 2*c*x^2))/(2*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) + (Sqrt[c]*(2*b - Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*
Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/2)*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*
(2*b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/2
)*Sqrt[b + Sqrt[b^2 - 4*a*c]])

________________________________________________________________________________________

Rubi [A]  time = 0.259512, antiderivative size = 221, normalized size of antiderivative = 1., number of steps used = 4, number of rules used = 3, integrand size = 18, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.167, Rules used = {1119, 1166, 205} \[ -\frac{x \left (b+2 c x^2\right )}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}+\frac{\sqrt{c} \left (2 b-\sqrt{b^2-4 a c}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (\sqrt{b^2-4 a c}+2 b\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{\sqrt{b^2-4 a c}+b}} \]

Antiderivative was successfully verified.

[In]

Int[x^2/(a + b*x^2 + c*x^4)^2,x]

[Out]

-(x*(b + 2*c*x^2))/(2*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) + (Sqrt[c]*(2*b - Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*
Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/2)*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*
(2*b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/2
)*Sqrt[b + Sqrt[b^2 - 4*a*c]])

Rule 1119

Int[((d_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[(d*(d*x)^(m - 1)*(b + 2*c*
x^2)*(a + b*x^2 + c*x^4)^(p + 1))/(2*(p + 1)*(b^2 - 4*a*c)), x] - Dist[d^2/(2*(p + 1)*(b^2 - 4*a*c)), Int[(d*x
)^(m - 2)*(b*(m - 1) + 2*c*(m + 4*p + 5)*x^2)*(a + b*x^2 + c*x^4)^(p + 1), x], x] /; FreeQ[{a, b, c, d}, x] &&
 NeQ[b^2 - 4*a*c, 0] && LtQ[p, -1] && GtQ[m, 1] && LeQ[m, 3] && IntegerQ[2*p] && (IntegerQ[p] || IntegerQ[m])

Rule 1166

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
 q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 205

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]*ArcTan[x/Rt[a/b, 2]])/a, x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rubi steps

\begin{align*} \int \frac{x^2}{\left (a+b x^2+c x^4\right )^2} \, dx &=-\frac{x \left (b+2 c x^2\right )}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}+\frac{\int \frac{b-2 c x^2}{a+b x^2+c x^4} \, dx}{2 \left (b^2-4 a c\right )}\\ &=-\frac{x \left (b+2 c x^2\right )}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\left (c \left (1+\frac{2 b}{\sqrt{b^2-4 a c}}\right )\right ) \int \frac{1}{\frac{b}{2}+\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx}{2 \left (b^2-4 a c\right )}+\frac{\left (c \left (2 b-\sqrt{b^2-4 a c}\right )\right ) \int \frac{1}{\frac{b}{2}-\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx}{2 \left (b^2-4 a c\right )^{3/2}}\\ &=-\frac{x \left (b+2 c x^2\right )}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}+\frac{\sqrt{c} \left (2 b-\sqrt{b^2-4 a c}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (2 b+\sqrt{b^2-4 a c}\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b+\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{b+\sqrt{b^2-4 a c}}}\\ \end{align*}

Mathematica [A]  time = 0.442591, size = 222, normalized size = 1. \[ \frac{-b x-2 c x^3}{2 \left (b^2-4 a c\right ) \left (a+b x^2+c x^4\right )}-\frac{\sqrt{c} \left (\sqrt{b^2-4 a c}-2 b\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{b-\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{b-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (\sqrt{b^2-4 a c}+2 b\right ) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{c} x}{\sqrt{\sqrt{b^2-4 a c}+b}}\right )}{\sqrt{2} \left (b^2-4 a c\right )^{3/2} \sqrt{\sqrt{b^2-4 a c}+b}} \]

Antiderivative was successfully verified.

[In]

Integrate[x^2/(a + b*x^2 + c*x^4)^2,x]

[Out]

(-(b*x) - 2*c*x^3)/(2*(b^2 - 4*a*c)*(a + b*x^2 + c*x^4)) - (Sqrt[c]*(-2*b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]
*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/2)*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]
*(2*b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*(b^2 - 4*a*c)^(3/
2)*Sqrt[b + Sqrt[b^2 - 4*a*c]])

________________________________________________________________________________________

Maple [A]  time = 0.232, size = 342, normalized size = 1.6 \begin{align*}{\frac{x}{8\,ac-2\,{b}^{2}} \left ({x}^{2}+{\frac{b}{2\,c}}-{\frac{1}{2\,c}\sqrt{-4\,ac+{b}^{2}}} \right ) ^{-1}}+{\frac{c\sqrt{2}b}{4\,ac-{b}^{2}}{\it Artanh} \left ({cx\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}-{\frac{c\sqrt{2}}{8\,ac-2\,{b}^{2}}{\it Artanh} \left ({cx\sqrt{2}{\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( -b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{x}{8\,ac-2\,{b}^{2}} \left ({x}^{2}+{\frac{1}{2\,c}\sqrt{-4\,ac+{b}^{2}}}+{\frac{b}{2\,c}} \right ) ^{-1}}+{\frac{c\sqrt{2}b}{4\,ac-{b}^{2}}\arctan \left ({cx\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{-4\,ac+{b}^{2}}}}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}}+{\frac{c\sqrt{2}}{8\,ac-2\,{b}^{2}}\arctan \left ({cx\sqrt{2}{\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \right ){\frac{1}{\sqrt{ \left ( b+\sqrt{-4\,ac+{b}^{2}} \right ) c}}}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^2/(c*x^4+b*x^2+a)^2,x)

[Out]

1/2/(4*a*c-b^2)*x/(x^2+1/2*b/c-1/2/c*(-4*a*c+b^2)^(1/2))+c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((-b+(-4*a*c
+b^2)^(1/2))*c)^(1/2)*arctanh(x*c*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b-1/2*c/(4*a*c-b^2)*2^(1/2)/((-b+
(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(x*c*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2))+1/2/(4*a*c-b^2)*x/(x^2+1/2
/c*(-4*a*c+b^2)^(1/2)+1/2*b/c)+c/(4*a*c-b^2)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arcta
n(x*c*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))*b+1/2*c/(4*a*c-b^2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*a
rctan(x*c*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2))

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2/(c*x^4+b*x^2+a)^2,x, algorithm="maxima")

[Out]

-1/2*(2*c*x^3 + b*x)/((b^2*c - 4*a*c^2)*x^4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2) - 1/2*integrate((2*c*x^2
- b)/(c*x^4 + b*x^2 + a), x)/(b^2 - 4*a*c)

________________________________________________________________________________________

Fricas [B]  time = 1.7401, size = 3623, normalized size = 16.39 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2/(c*x^4+b*x^2+a)^2,x, algorithm="fricas")

[Out]

-1/4*(4*c*x^3 + sqrt(1/2)*((b^2*c - 4*a*c^2)*x^4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2)*sqrt(-(b^3 + 12*a*b*
c + (a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5
*c^3))/(a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))*log((3*b^2*c + 4*a*c^2)*x + 1/2*sqrt(1/2)*(b^5 -
8*a*b^3*c + 16*a^2*b*c^2 - (a*b^8 - 8*a^2*b^6*c + 128*a^4*b^2*c^3 - 256*a^5*c^4)/sqrt(a^2*b^6 - 12*a^3*b^4*c +
 48*a^4*b^2*c^2 - 64*a^5*c^3))*sqrt(-(b^3 + 12*a*b*c + (a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sq
rt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)
)) - sqrt(1/2)*((b^2*c - 4*a*c^2)*x^4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2)*sqrt(-(b^3 + 12*a*b*c + (a*b^6
- 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b
^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))*log((3*b^2*c + 4*a*c^2)*x - 1/2*sqrt(1/2)*(b^5 - 8*a*b^3*c +
 16*a^2*b*c^2 - (a*b^8 - 8*a^2*b^6*c + 128*a^4*b^2*c^3 - 256*a^5*c^4)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2
*c^2 - 64*a^5*c^3))*sqrt(-(b^3 + 12*a*b*c + (a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6
- 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))) + sqrt(1
/2)*((b^2*c - 4*a*c^2)*x^4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2)*sqrt(-(b^3 + 12*a*b*c - (a*b^6 - 12*a^2*b^
4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2
*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))*log((3*b^2*c + 4*a*c^2)*x + 1/2*sqrt(1/2)*(b^5 - 8*a*b^3*c + 16*a^2*b*c
^2 + (a*b^8 - 8*a^2*b^6*c + 128*a^4*b^2*c^3 - 256*a^5*c^4)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a
^5*c^3))*sqrt(-(b^3 + 12*a*b*c - (a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^
4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))) - sqrt(1/2)*((b^2*c
 - 4*a*c^2)*x^4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2)*sqrt(-(b^3 + 12*a*b*c - (a*b^6 - 12*a^2*b^4*c + 48*a^
3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2*b^4*c + 48
*a^3*b^2*c^2 - 64*a^4*c^3))*log((3*b^2*c + 4*a*c^2)*x - 1/2*sqrt(1/2)*(b^5 - 8*a*b^3*c + 16*a^2*b*c^2 + (a*b^8
 - 8*a^2*b^6*c + 128*a^4*b^2*c^3 - 256*a^5*c^4)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^4*b^2*c^2 - 64*a^5*c^3))*sq
rt(-(b^3 + 12*a*b*c - (a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3)/sqrt(a^2*b^6 - 12*a^3*b^4*c + 48*a^
4*b^2*c^2 - 64*a^5*c^3))/(a*b^6 - 12*a^2*b^4*c + 48*a^3*b^2*c^2 - 64*a^4*c^3))) + 2*b*x)/((b^2*c - 4*a*c^2)*x^
4 + a*b^2 - 4*a^2*c + (b^3 - 4*a*b*c)*x^2)

________________________________________________________________________________________

Sympy [A]  time = 3.74574, size = 298, normalized size = 1.35 \begin{align*} \frac{b x + 2 c x^{3}}{8 a^{2} c - 2 a b^{2} + x^{4} \left (8 a c^{2} - 2 b^{2} c\right ) + x^{2} \left (8 a b c - 2 b^{3}\right )} + \operatorname{RootSum}{\left (t^{4} \left (1048576 a^{7} c^{6} - 1572864 a^{6} b^{2} c^{5} + 983040 a^{5} b^{4} c^{4} - 327680 a^{4} b^{6} c^{3} + 61440 a^{3} b^{8} c^{2} - 6144 a^{2} b^{10} c + 256 a b^{12}\right ) + t^{2} \left (- 12288 a^{4} b c^{4} + 8192 a^{3} b^{3} c^{3} - 1536 a^{2} b^{5} c^{2} + 16 b^{9}\right ) + 16 a^{2} c^{3} + 24 a b^{2} c^{2} + 9 b^{4} c, \left ( t \mapsto t \log{\left (x + \frac{16384 t^{3} a^{5} c^{4} - 8192 t^{3} a^{4} b^{2} c^{3} + 512 t^{3} a^{2} b^{6} c - 64 t^{3} a b^{8} - 128 t a^{2} b c^{2} - 16 t a b^{3} c - 4 t b^{5}}{4 a c^{2} + 3 b^{2} c} \right )} \right )\right )} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**2/(c*x**4+b*x**2+a)**2,x)

[Out]

(b*x + 2*c*x**3)/(8*a**2*c - 2*a*b**2 + x**4*(8*a*c**2 - 2*b**2*c) + x**2*(8*a*b*c - 2*b**3)) + RootSum(_t**4*
(1048576*a**7*c**6 - 1572864*a**6*b**2*c**5 + 983040*a**5*b**4*c**4 - 327680*a**4*b**6*c**3 + 61440*a**3*b**8*
c**2 - 6144*a**2*b**10*c + 256*a*b**12) + _t**2*(-12288*a**4*b*c**4 + 8192*a**3*b**3*c**3 - 1536*a**2*b**5*c**
2 + 16*b**9) + 16*a**2*c**3 + 24*a*b**2*c**2 + 9*b**4*c, Lambda(_t, _t*log(x + (16384*_t**3*a**5*c**4 - 8192*_
t**3*a**4*b**2*c**3 + 512*_t**3*a**2*b**6*c - 64*_t**3*a*b**8 - 128*_t*a**2*b*c**2 - 16*_t*a*b**3*c - 4*_t*b**
5)/(4*a*c**2 + 3*b**2*c))))

________________________________________________________________________________________

Giac [F(-2)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Exception raised: NotImplementedError} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2/(c*x^4+b*x^2+a)^2,x, algorithm="giac")

[Out]

Exception raised: NotImplementedError